FRW cosmology in Milgrom's bimetric theory of gravity 



Timothy Clifton 1 * and Thomas G. Zlosnik 2 t 
1 Department of Astrophysics, University of Oxford, UK. 
2 Perimeter Institute of Theoretical Physics, Waterloo, Canada. 

(Dated: May 21, 2010) 

We consider spatially homogeneous and isotropic Friedmann-Robertson- Walker (FRW) solutions 
of Milgrom's recently proposed class of bimetric theories of gravity. These theories have two different 
regimes, corresponding to high and low acceleration. We find simple power-law matter dominated 
solutions in both, as well as solutions with spatial curvature, and exponentially expanding solutions. 
In the high acceleration limit these solutions behave like the FRW solutions of General Relativity, 
with a cosmological constant term that is of the correct order of magnitude to explain the observed 
accelerating expansion of the Universe. We find that solutions that remain in the high acceleration 
regime for their entire history, however, require non-baryonic dark matter fields, or extra interaction 
terms in their gravitational Lagrangian, in order to be observationally viable. The low acceleration 
regime also provides some scope to account for this deficit, with solutions that differ considerably 
from their general relativistic counterparts. 
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I. INTRODUCTION 

A natural way to extend General Relativity is via the 
introduction of a second dynamical rank 2 metric tensor 
field, g a b. The so called 'bimetric' theories that result 
have a rich phenomenology (see, for example, [1-4]), and 
a long history 5-7]. They represent an application to 
gravity of the Yang-Mills approach to gauge theories, and 
allow for new and interesting behaviour. 

A new class of these theories has recently been intro- 
duced by Milgrom @, Q , in which the interaction term 
between the two metric fields is constructed from a ten- 
sor defined as the difference of the connections associated 
with each of them. The stated purpose of this approach 
is to produce a weak-field limit of the form 

V 2 *-V- [/(Z|V*AH)Vtfjv] =4rrGp, (1) 

where &n obeys V 2 ^n = kxGp, and where time-like 
test particles obey a force law x ~ — V^F. Here \& is a 
gravitational potential, f ix) is a function to be specified, 
p is the energy density and I is a constant with units of 
length. The existence of gravitational fields of this type 
would be of considerable interest for astrophysics [l(J □J, 
and the goal of producing a viable relativistic formulation 
for theories of this kind has been pursued for some time 

Let us now try to understand why a relativistic formu- 
lation of ([T]) is so important. If we associate an energy- 
momentum tensor, T a b, with the right-hand side of (TTJ), 
then 

p = T ab n a n b , (2) 
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where n a is the unit vector normal to surfaces of constant 
t. Energy-momentum conservation then tells us that test 
particles follow geodesies of the metric with respect to 
which T a b is covariantly conserved. For a perturbed met- 
ric of the form 

ds 2 = -(l + 2*)dt 2 + (l-2<f>)6 ij dx i dx j (3) 

this corresponds to x ~ — V^, as required. In order to 
gain a set of field equations that are invariant with re- 
spect to coordinate transformations, and that reduce to 
(HJ) for the metric ([3|), we must now construct a tensor 
quantity that can be associated with the left-hand side of 
(fTj). We immediately encounter a problem though: The 
quantity |V V I'| cannot be formed from any local curva- 
ture invariants of (j3|) alone. One may alternatively try 
and recover (TTJ from nonlocal functions of the space-time 
curvature [2lj T However, such attempts have prompted a 
no-go theorem [22j regarding the phenomenological via- 
bility of all theories that are constructed from the space- 
time geometry and matter fields only. New structure is 
therefore required, beyond a single space-time metric. 

In the absence of a relativistic theory that is well moti- 
vated by other considerations, and that reduces to ([TJ in 
the appropriate limit, it is not clear what this additional 
structure should be. We are then left without any way 
of knowing when JTJ) should be considered valid, or when 
we should prefer, for example, a Newtonian description. 
We can, however, use observations of relativistic effects 
beyond Newtonian order to gain insight. In this regard, 
cosmological solutions provide us with a very useful probe 
as they express the full non-linearity of the relativistic 
field equations. Motivated by this, as well as the desire 
to understand further the cosmological consequences of 
bimetric theories of gravity in general, we consider here 
the FRW solutions of Milgrom's class of bimetric theo- 
ries. 

Of course, for any theory of gravity to be considered 
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viable it is now the case that it should be able to repro- 
duce all of the major probes of observational cosmology. 
These include the anisotropics in the cosmic microwave 
background, the matter power spectrum on large scales, 
Hubble diagrams that extend out to z ~ 1 , as well as the 
ratios of light elements from primordial nucleosynthesis. 
All of these require a detailed knowledge of FRW cosmol- 
ogy (as well as of cosmological perturbation theory, in the 
case of the first two). It is not our goal here to perform 
an exhaustive study of all of these areas, but rather to 
make a solid first step in understanding the underlying 
FRW solutions, and what they imply. 

The article proceeds as follows. In Section [TT| we de- 
scribe the theory and its field equations. The general 
form of the theory allows considerable flexibility, and so 
in Section [Til] we discuss and motivate the specific ver- 
sions of the theory that we will consider. In Section lTVl wo 
find the field equations in space-times with FRW symme- 
tries, and discuss their solutions in Section [V] In Section 
IVII we discuss the degree to which current probes of cos- 
mology can constrain these models, and in Section IVIII 
we conclude. 



II. FIELD EQUATIONS 

The theory we are considering has two metrics, g a b and 
(jab- The action for the gravitational part of Milgrom's 
theory can be written in terms of these two tensor fields 
as 



I = Ig+Ig+ If, 



where 



(4) 

(5) 
(6) 



are the Einstein- Hilbert terms for g a b and g a b, and Ii n t 
is an interaction term between them. In the equations 
above R is the Ricci scalar constructed from g a b, and R 
is the Ricci scalar constructed from g a b- 

In Milgrom's theory the interaction term, J mt , is speci- 
fied as a function of scalars formed from the rank-3 tensor 




be — 1 be ~ 1 6c> 



(7) 



where T a bc and f a bc are the metric connections of g a b 
and g a b, respectively, such that 



9ab:c — 0, 



(8) 
(9) 



where ; and : indicate covariant derivatives constructed 
from r a bc and T a bc - The general form of scalars that are 
quadratic in C a bc can then be written as 



i2 n beef n a 

1 ^ad ° ic u e/i 



(10) 



where Q ad ce ^ is built from g a b and g a b, and I is a constant 
with dimensions of length. In terms of these quantities 
we can then write 1 



-gaM(Q)d 4 x, 



(11) 



where M.{Q) and a = <j{k) are functions to be specified, 
and k — {g/g) 1 ^ 4 - 

In Appendix A we find the field equations for g a b can 
be written as 



PG a b + S a b = SnGTab 



(12) 



where 



S ab = 2 aM Q U (ab) + J (a b) - J c al 



SQ 



ad HCef r g C d g(l + n) M 
j-aT— ° he c ef p Mg a . 



and G a b and T a b are the Einstein tensor and the energy- 
momentum tensor, defined with respect to g a b in the 
usual way. We have also used the notation Mq = 
dM/dQ and J/ f = Q ad bcef C a bc , and taken a = n 2n . 
Similarly, the field equations for g a b are found to be 



aGab + S a b = SnGTab 



(13) 



where 
S a b — 



ct"" 1 Mt 



(ef); 



j d cf {r x ) de 



dcagfb 



zcef 



n + l 

a " n 



-Mg a b, 



-^M Q 5 ^ b - C y zc c% + /2 



and where G a b and T a b are the Einstein and energy- 
momentum tensors defined with respect to g a b- Here, 
and throughout, we write the inverse of g a b as (g~ 1 ) ab 
(defined such that (<7 _1 ) ac 3b c = S a b ). Any raising or low- 
ering of indices is otherwise always done with g a b- 

It should be noted that in (fT"2]) and (p~3|) above we have 
taken the matter fields described by T a b and T a b to be 
coupled to g a b and g a b, respectively 2 . We presume that 
we (as observers) are made from matter coupled to only 
one of these metrics, which we take to be the former with- 
out loss of generality. We further presume that there is 
no interaction between T a b and T a b, so that all observa- 
tions we make will be of the other matter fields coupled 
to g a b- Cosmological probes of the expanding Universe 
then give us direct information about the geometry of g a b 
only. 



1 Note that in terms of the formalism used by Milgrom a = /(k)/k, 
and Q = i 2 S. 

2 Although one could also concievably couple matter to a combi- 
nation of these metrics, we prefer to restrict ourselves here to the 
case considered by Milgrom in 
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In the equations derived above we have not explic- 
itly included any cosmological constant terms. However, 
from (fT2jl and (fl3|) it is clear that A is dynamically equiv- 
alent to a perfect fluid with p = —p. We therefore ac- 
count for any possible cosmological constants by allowing 
for them in T a b and T a b- 

III. SPECIFICATION OF THE THEORY 

In Section |n] we presented the field equations for the 
theory derived from the action specified by ((5]), dU and 
These equations provide constraints on, and specify 
the evolution of, the two dynamical rank two tensor fields 
g a b and g a b, but they also contain considerable freedom: 
The tensor Q ab cde and the function M (Q) have yet to be 
specified. In order to make progress in understanding the 
cosmological solutions of these theories, we will therefore 
now restrict ourselves to specific cases. 

Firstly, we will only consider Q a h Cd ^ — Qha 
that are formed from the metric tensor g a b- In this case 
the most general Q ab cde ^ that can be constructed is 

+ Ci S c a 5e b g d f + c 5 g a bg ce g df , (14) 

where c\, C2, C3, C4 and C5 are constants, and indices are 
assumed to be symmetrized appropriately. In this nota- 
tion, the 'concrete simple theory' of Milgrom is specified 
by C\ = 1, C2 = — 1 and C3 = C4 = C5 = 0. The form of 
Q ab cde ^ given in (|T4l) can then be seen to correspond to 
a generalization of this theory. 

Secondly, we will consider the function Ai(Q) to have 
the form specified by Milgrom in the low and high accel- 
eration limits of the theory. That is, when \Q\ > 1 we 
take 

M Q ~0, (15) 

and when |Q| <C 1 we take 

Mq^U- 1 ^. (16) 

The first of these corresponds to Newtonian gravitation 
in the non-relativistic limit (when c — > 00), and the sec- 
ond corresponds to modified gravitational dynamics in 



the low acceleration regime of the non-relativistic limit 
(when c — > 00 and I — > 0). In particular, when a + /3 = 
and T a b = one recovers (Q]) in the weak field limit, with 
/ = M.q [|[. The weak field limit with more general a 
and (3, and with T a b, has been explored in [§]. 

Here we will not be concerned, for the most part, with 
the transitional behaviour between these regimes. It is 
not clearly specified by the weak field limit of the theory, 
and is presumed to be highly sensitive to the particular 
form of A4(Q) that is chosen. 

IV. FRW COSMOLOGY 

We will now consider the cosmological solutions of the 
theory discussed above. Imposing FRW symmetries on 
g a b and g a b leads to the line-elements 

ds 2 = g ab dx a dx b = -a(r) 2 rfr 2 + a(r) 2 h^dx 1 dx\ (17) 

and 

ds 2 = g a bdx a dx b = -X{t) 2 <It 2 + Y{T) 2 h ll dx l dx j , (18) 

where hij and hij are the metrics of Euclidean 3-spaces 
of constant curvature, k and k, respectively. 

The reason for introducing the function X(t) in the 
T-r component of g a b is that we now only have one co- 
ordinate freedom of the form r — > /(r), but two t-t 
components, in the two different metrics. It is therefore 
not possible to redefine r to absorb X(t) without intro- 
ducing a second function into the t-t component of g a b- 
The most general expression of FRW geometry for both 
metrics must therefore contain three functions of t: <x(t), 
X(t) and Y(t). 

Using the geometry specified by (jTTJ) and (|18[) we can 
now calculate the quantities that appear in the field equa- 
tions (TT^]) and (fTB"]) . Of particular interest is the con- 
straint equation, which contains only first derivatives of 
a, X and Y . This equation is the analogue of the Fried- 
mann equation in general relativistic FRW cosmology, 
and is derived in Appendix B. As an example, for Mil- 
grom's 'concrete simple theory' with k = k (which we 
will motivate later) we find 
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3Z 2 
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and where p and p are the energy densities of the perfect 
fluids associated with the energy-momentum tensors T ao 
and Tab, respectively. Over-dots denote differentiation 
with respect to r. More general expressions for the con- 
straint equation, and for the other quantities involved in 
the field equations, are given in Appendix C. 

The field equations (|12p and (TT3")) also provide second- 
order evolution equations for the three variables a, X and 
Y. However, these expressions are very lengthy, and so 
we choose not to reproduce them explicitly here. 

As usual, the energy-momentum tensors obey conser- 
vation equations. For non-interacting fluids these equa- 
tions read 



where 



and 



p + 3-{p + p) = 



(19) 



(20) 



where p and p are the pressures of the two ideal fluids 3 . 
It is convenient to treat these fluids as being polytropic, 
with constant equations of state w and w defined by p = 
wp and p — wp. Equations (fT9"|) and (|2U1) then give p cx 

-3(i+t«) ^ p y-3(i+«i) 



V. SOLUTIONS 

We will now look for solutions to the field equations in 
both the 'high acceleration' limit, when \Q\ 1, and the 
'low acceleration' limit, when \Q\ -C 1. 



A. High acceleration limit 

For |Q| 3> 1 we have Mq ~ 0, so the constraint and 
evolution equations for g ao are given by 

a 2 _ 8nGp k Mo(l + n)a 



a a 

n3 ^4 



3/3 a 2 3/31 2 
4irG(p + 3p) M a (l + n)a 



3/3 



3/8/ 2 



. (22) 



while the constraint and evolution equations for g a b are 
given by 



Y 2 



X 2 Y 2 
Y XY 



8-irGp k M. na « 
3a ~ Y 2 + 



3a/ 2 



(23) 



X 2 Y X 3 Y 



AnG(p + 'dp) + Mpna^ ^ 



3o 



3a/ 2 



3 It is, in principle, possible to have an interaction between the two 
fluids such that energy-momentum could be exchanged between 
them. We do not consider this possibility here, but rather treat 
the two fluids as being non-interacting, as is usual in cosmology. 



\f\ — kr 2 a 4 

7f _ kr 2 xy^ 



(25) 



and Mo is the constant value of A4(x) when i» 1. It 
is now clear that one or more of the following conditions 
must be true: 

(i) The value of Mq is 0. 

(ii) The value of n is or — 1. 

(Hi) The value of k is equal to that of k. 

If none of these conditions are met then the terms con- 
taining a have an r dependence that cannot be accom- 
modated by any of the other terms. The first two of these 
are conditions on the theory. If either, or both, of these 
conditions are met then there can exist FRW solutions in 
the two metrics which have different spatial curvatures. 
If neither (i) nor (ii) are satisfied, then the only FRW 
solutions that exist are those in which the spatial curva- 
tures of the two FRW metrics are the same. 

If condition (i) is met then the Friedmann equations 
above can be seen to be unchanged from their general 
relativistic form. The two scale factors a and Y are 
completely decoupled, and obey constraint and evolution 
equations that are exactly the same as an FRW geometry 
in General Relativity. 

If condition (ii) is true (but (i) is not) then either a or 
Y is driven by an additional term in its field equations 
that acts in exactly the same as a cosmological constant. 
For n — this term can be seen to appear in the Fried- 
mann equations for a, which are otherwise identical to 
their general relativistic counterparts. The additional ef- 
fective cosmological constant is given by 



A 



eff 



Mo 

' f3l 2 ■ 



(26) 



If this is the case then Y obeys field equations that are 
unchanged from a general relativistic FRW solution. If 
n = — 1 then the situation is reversed, with Y being 
driven by an additional term, while the Friedmann equa- 
tions for a are unchanged. 

If the theory does not satisfy conditions (i) or (ii) then 
the only FRW solutions that exist have k — k. The 
value of a is then independent of r, and a obeys a set 
of Friedmann equations with an additional effective fluid 
whose energy density and pressure are given by 



Peff 



-Peff 



M (l + n) ( a 4 



8ttG/ 2 \XY 3 



(27) 



This effective fluid has an equation of state w e f / = — 1 , 
and so we must have p e // =constant, and therefore 
a oc A 1 / 4 !" 3 / 4 . In this case the other set of Friedmann 
equations also has an additional effective fluid, with 



Peff 



-Peff 



Mon ( a 4 



8irGl 2 \ XY 3 



(28) 
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We therefore have that one of the scale factors is driven 
by an additional positive effective cosmological constant, 
while the other is driven by a negative one (unless < 
n < — 1, in which case both effective cosmological con- 
stant terms can have the same sign). 

It has been shown above that in all possible cases 
the two FRW geometries effectively decouple from each 
other. Their evolution is then specified by a set of equa- 
tions that are identical to the Friedmann equation of 
General Relativity, but with the possible addition of an 
effective cosmological constant term. Further more, un- 
less condition (i) or (ii) is satisfied, we must also have 
k = k and a cx X^Y 3 ^. 



Low acceleration limit 



where p and q are given by 



I - 3w 



P 



(33) 



(I + 2w + 8nw — 3w — 8nw — 6ww) 

_ l-3w , 34 > 

' (I + 2w + 8nw — 3w — 8nw — 6ww) ' 

and where w and w are the equations of state of the two 
fluids, as defined in Section [IVJ These solutions are sig- 
nificantly different from the corresponding solutions in 
General Relativity, and some examples, for cosmologi- 
cally interesting equations of state, are given in Table 
HI In this table, and throughout, we use the time co- 
ordinates t and t to correspond to the proper time of 
comoving observers in each geometry, so that t = J adr 
and t = J Xdr. 



The situation in the low acceleration limit is more com- 
plicated. In this case we have |Q| <C 1 and 



M ~ sign(Q) 



3/4 



Mi 



(29) 



-1/4 



Here Mi is an unspecified con- 



so that Mq ~ 
stant. 

Firstly, let us consider the situation with k = k = 
and Mi — (we will return to the situation with non- 
zero spatial curvature, and Mi ^ 0, below). In this case 
it can seen from Appendix C that we automatically have 
Q = Q( T ) and cr = o"(t)i an d that all dependence on r 
vanishes from the field equations. Now, just as one can 
take c — > oo to find the non-relativistic limit of General 
Relativity, here we take I — > to find the low acceleration 
limit of the present theory. Such a limit suppresses the 
contribution of G a b and G a b to the left hand side of the 
field equations, relative to S a b and S a b, so that (|I2p and 
(fT3)) become 



and 



S a b — 8irGT a b 



S a b — 8nGT a b 



(30) 



(31) 



These equations are second-order, and provide a system 
of constraint and evolution equations. Assuming a per- 
fect fluid form for both energy-momentum tensors we 
then find the following matter dominated power-law so- 
lutions 4 : 



and 



X oc Y oc r q 



(32) 



4 By 'matter dominated' we mean that the left-hand side of the 
field equations scales in the same way as the right-hand side. 
This is in contrast to 'vacuum dominated', in which the left- 
hand side scales independently of the right-hand side, such that 
the influence of the matter fields on the space-time dynamics is 
negligible. 



Fluids 



dust & dust 
dust & rad. 
rad. & dust 
rad. & rad. 
scalar & scalar 
A & A 



a(t) 



.1/2 



constant 



3+SnC-3C 
t l/4 

exponential 



Y(t) 



1/2 



£(3-8t.) 

constant 
t° 

£1/4 

exponential 



GR 



t 2/3 

t l/2 
t l/3 

exponential 



TABLE I: The functional form of matter dominated power- 
law solutions in the low acceleration limit of Milgrom's bimet- 
ric theory and General Relativity (GR). The fluids considered 
are dust (w — 0), radiation (w = 1/3), scalar fields (w — 1) 
and A (w = — 1). In the left hand column, the first fluid is 
coupled to g a b, and the second fluid to g a b- In the radiation 
& radiation case, the solutions are not specified uniquely, and 
hence there is an arbitrary constant C involved. 

It can be seen that the functional form of these solu- 
tions is, in general, not only dependent on the fluid that 
is coupled to the metric in question, but also to the fluid 
coupled to the other metric, as well as the theory, via 
the parameter n. This mutual dependence should be ex- 
pected as the two metrics are non-minimally coupled to 
each other in the action. The evolution of the two met- 
rics can also be asymmetric under the exchange of the 
two fluids. This can be seen from the second and third 
rows of Table HI Again, this is expected, as Q (as speci- 
fied in Section IIII[) is asymmetric under the exchange of 
g a b and g a b- Once these power-law solutions have been 
assumed, the field equations then reduce to a set of al- 
gebraic relations between the Cj's, n, a, /3, I, p(to) and 
p(to), as well as the constants of proportionality in (|32l) . 
where to is the present age of the Universe. 

It is interesting to note that in some cases it is possible 
for one of the metrics to be static. This is exemplified in 
Table U by the two radiation & dust solutions. In fact, it 
can be seen from ([33]) and (f34|) that if the matter coupled 
to the second metric is radiation, then the first metric 
will be static (unless the denominator diverges, as is the 
case with radiation & radiation). It is also interesting 
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to note that in the radiation & radiation case the two 
scale factors are not uniquely determined, and can only 
be specified up to an arbitrary constant. Clearly the 
behaviour of these solutions is quite different to the usual 
general relativistic FRW solutions, which are shown in 
the fourth column of Table U for comparison. 

A special case that should be noted is when X equals 
Y , rather than simply being proportional to it. In this 
case the solutions that are obtained are quite different to 
the general case, discussed above. What happens is that 
the contributions from any terms that involve Q vanish. 
The contribution of G a b and G a b can then no longer be 
considered to be suppressed, as the terms they were sup- 
pressed with respect to no longer exist. The system of 
equations then reduces to exactly those considered pre- 
viously, in the high acceleration regime. This is, in fact, 
the case that was considered by Milgrom in Q. 

Now let us consider non-zero spatial curvature. It can 
be seen from Appendix C that a and Q are, in general, 
functions of both r and r, when k and k take arbitrary 
values. Such a dependence on r is problematic, as the free 
functions in the line-elements (IT71) and (jT5J) are all func- 
tions of r only. This is a generalization of the problem 
involving the r dependence of a found in the high ac- 
celeration regime, and although the dependence is more 
complicated, the solution is the same: a and Q reduce 
to functions of r if k = k. It can then be seen that not 
only does all r dependence drop out of a, Q, S a b and 
S a b, but that these functions also become independent of 
k and k (up to possible overall multiplicative factors of 
(1 — kr 2 )). The only place that non-zero k has any effect 
in the field equations is therefore through its contribu- 
tion to G a b and G a b, where it appears as an additional 
constant. 

Now, although the contribution of the fc-independent 
terms in G a b are suppressed by a power of I, and can 
therefore be safely ignored in the low acceleration regime, 
the same cannot be said for k. This term is a constant, 
and does not have the potential to become arbitrarily 
small at late-times, in contrast to the other terms in G a b 
and S a b- We therefore cannot ignore it, as even if it is 
negligibly small at early times, due to suppression by a 
power of I, it can still become influential at late-times. 
However, as the influence of k and k is limited to its ap- 
pearance in the Einstein tensor, when k = k, its influence 
can be straightforwardly accounted for: It acts in exactly 
the same way as an additional fluid with w = —1/3 cou- 
pled to g a b, and a fluid with w = —1/3 coupled to g ao . We 
can therefore calculate its influence at late-times using 
equations (|3"3")l and (|34p. above. We find that the corre- 
sponding power-law solution has <zcxAcxY"cxt 3 cx i 3 / 4 , 
independent of n. This should be contrasted with the 
usual a oc t in General Relativity. 

Next, let us consider the effect of non-zero M.\. We 
find that the inclusion of such a term requires a cx t _1 , 
if power-law solutions are to exist, which correspond to 
exponential evolution in the proper time of comoving ob- 
servers, t. The solutions that exist for g ao then depend 



on the matter content of the theory. If we couple a fluid 
with equation of state w to g a b then we find the power- 
law solution 



X cx Y cx r" 



OC f 1+ 3(i+™> , 



(35) 



while if we couple a fluid with equation of state w to g ao 
we find 



+ A 4(l+n) 

X txY (xt ci+4«-!wo cx i3<i+d) . 



(36) 



Power-law solutions with fluids coupled to both metrics 
do not, in general, exist, unless they result in a cx r _1 
anyway (such as, for example, two A terms). 

Finally, let us consider what happens on approach to 
a vacuum. For monotonically expanding geometries flu- 
ids with higher equations of state give way to those with 
lower equations of state, such that dust domination fol- 
lows radiation domination, and so forth, in the usual way. 
If a fluid with equation of state w = — 1 (or less) exists 
then the total energy density will eventually reach a con- 
stant (or start to increase), as all other fluids become 
negligible in comparison. If no such fluid exists, then the 
total energy density will decrease monotonically forever. 
If this happens then the right-hand side of the field equa- 
tions can eventually drop below the previously neglected 
terms in the Einstein tensor. For p(l + 3u>) > 2 this 
will mean that the matter becomes subdominant, and the 
cosmological evolution will be vacuum dominated (deter- 
mined by the dynamics of the vacuum alone) . In this case 
the only power-law solutions that exist go like t , which 
corresponds to exponential expansion. For p(l + 3u>) < 2 
the right-hand side of the field equations will always be 
dominant over the neglected terms in the Einstein ten- 
sor, and matter dominated power-law evolution, as de- 
scribed above, can continue indefinitely into the future. 
It should be noted that if spatial curvature is non-zero, 
and acts like a fluid with w = —1/3, then the condition 
p(l + 3w) < 2 is always satisfied. A curvature dominated 
power-law solution, as found above, can therefore always 
last for an indefinitely long period, provided there are no 
fluids with equations of state w < —1/3. 



VI. VIABILITY AS A MODEL OF THE 
UNIVERSE 

We will now consider the viability of the cosmological 
solutions found above as models of the observable Uni- 
verse. Our criterion for viability will be whether or not 
it is possible to account for the major probes of obser- 
vational cosmology, which we take to be the primordial 
synthesis of light elements, the position of the first acous- 
tic peak in the CMB angular power spectrum, the growth 
of structure, and the late-time accelerating expansion of 
the Universe. For later convenience, let us now define the 
density fraction of a fluid i to be 



2 ' 



(37) 
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where pi is the energy density of the fluid on a homo- 
geneous hypersurface of age t , and H is the locally 
measured value of the Hubble parameter (in whichever 
geometry the observers are taken to be coupled to). 

First consider the primordial synthesis of light ele- 
ments. This process requires a cosmological evolution 
that is close to the standard form [23[, a ~ t 1 / 2 , and is 
highly sensitive to the ratio of photons to baryons, r\. In 
particular, at the end of nucleosynthesis (t ~ 200s) it can 
be shown that the observed abundances of light elements 
imply rj = (5.5 ± 0.5) x 10~ 10 [13] ■ If we now combine 
this with a CMB photon temperature of 2.725 ± 0.001K, 
and assume that the expansion of the Universe has been 
adiabatic, then we are left with the constraint 



= 0.020 ±0.002, 



(38) 



to 95% confidence, where h = iJoAlOOkms^Mpc" 1 ). 
This is a tight constraint on the matter content of the 
Universe, and one that we expect to be applicable to any 
viable model. 

Next let us consider the positions of the acoustic peaks 
of the CMB. In general relativistic FRW cosmology their 
angular extent on the sky, 9, is essentially determined by 
two factors: The acoustic horizon at last scattering, Ri s , 
and the angular diameter distance to the last scattering 
surface, d* A , such that 



Ris 



(39) 



Here the situation may be somewhat different. The for- 
mation of linear structure up to last scattering can be 
modified, as the weak-field limit can now be altered from 
its usual Newtonian form. This may affect the scale of 
correlations on the last scattering surface, but the pro- 
jection of that surface onto our sky will be unaffected by 
this, and will therefore be sensitive to the intervening ge- 
ometry of space-time in the usual way. We should there- 
fore expect an equation similar in form to (|39l) , but with 
Ris now replaced by i?; s , referring to the new quantity 
that sets the scale of correlations on the last scattering 
surface (no longer necessarily the acoustic horizon). 

A reliable calculation of Ri s will require a detailed anal- 
ysis of linear perturbation theory, and pre-recombination 
physics, within the frame- work of the current theory. For 
our current purposes, we will simply assume Ri s has been 
calculated. Using the well known result that spatial cur- 
vature produces a shift in the angular scale of the first 
acoustic peak according to (1 — fifc) -9 / 20 , and given the 
compatibility of the observed CMB with a spatially flat 
general relativistic FRW cosmology, we can then specu- 
late that, in the context of the present theory, the posi- 
tion of the first acoustic peak is likely to imply an effective 
spatial curvature of 



fin 
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This can be seen to reduce to ftk — when Ri s ~ R. 

Beyond the sensitivity of CMB observations to the cur- 
vature of the Universe, the growth of structure also de- 
pends on the evolution of the Universe. This is true for 
both the linear structure that we observe at last scat- 
tering, as well as the subsequent growth of non-linear 
structures such as galaxies, and clusters of galaxies. In 
most cosmological models, it is usually required that one 
should have a period of cosmological expansion where 
a t 2 / 3 in order for structure to form. Here the situation 
may be somewhat different, as, once again, the weak-field 
limit is modified from its usual Newtonian form. Never 
the less, we will consider it preferable if the cosmological 
solutions we have found can be shown to exhibit a period 
where a or X ~ t 2 l z . 

Finally, let us consider supernova observations. The 
Hubble diagrams constructed from these events allow us 
to determine the late-time evolution of the Universe. In 
the standard cosmological model, with the assumption of 
spatial flatness, they provide strong, direct evidence for a 
cosmological constant, or vacuum energy, with Q\ ~ 0.7 
[25ll26|. If one permits non-zero spatial curvature, as will 
be required in these models if Ri s ^ Ri Sl then the result- 
ing bounds on Q A change. However, it can still be shown 
that there exists good evidence for a non-zero cosmo- 
logical constant, independent of the value of fife. Using 
bounds derived from the 414 supernovae of the "Union" 
dataset one can make the approximate statement that 
27 1 



0.3 < n A < 1.4, 



(41) 



to 3cr, independent of the amount of spatial curvature 5 . 
Any viable model of the Universe should therefore include 
a late-time period of accelerating expansion. 

Having discussed the constraints available from obser- 
vational cosmology, let us now consider their applicability 
to the solutions found in Section IVl 



A. High acceleration limit 

It was shown above that when \Q\ ^> 1 the gravita- 
tional field equations reduce to their usual general rel- 
ativistic form, with a cosmological constant term. For 
a cosmological solution with dust and radiation we will 
therefore have periods of cosmological evolution with 
a ~ i 1 / 2 and i 2 / 3 , when these two fluids dominate the cos- 
mological dynamics, respectively. This is in good keeping 
with the requirements of primordial nucleosynthesis and 
structure formation. 

If we apply the primordial nucleosynthesis bound (|38[) 
on Qb, and take this to make up the entire contribu- 



(40) 



The precise bounds achievable depend on exactly which data one 
chooses to include, and how it is treated. 
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tion to the dust content of the Universe, then the super- 
novae results [27jj imply a low value of fl\ ~ 1/3, that is 
only compatible with observations at around the 3cr level. 
What is more, without any other matter content this 
would imply ttk ~ 2/3, which would only be compatible 
with the position of the first acoustic peak of the CMB 
angular power spectrum if Ri s ~ 0.6i?; s . Such a large 
change of scale seems unlikely from previous studies us- 
ing similar theories (2|| [29[ . We therefore find that some 
additional non-baryonic dust-like degrees of freedom are 
necessary, either in the form of cosmological dark matter, 
or from some other sector of the theory. 

An effective cosmological constant term, Mo, is also 
present in (|2ip and (|22j) . This term has the attractive 
feature of being able to reproduce the value of the cos- 
mological constant observed in supernova data with O(l) 
values of Aio and n, for astrophysically interesting values 
of I. It also provides a lower limit for \Q\, as Q — >constant 
in the limit that Mq dominates. This constant is zero 
in the special case X = Y, and, as already discussed, in 
this case the dynamics in the low acceleration limit be- 
come identical to those of the high acceleration regime. 
More generally we find that in the limit a —> aor" 1 , 
X -> Xqt^ 1 and Y -» Yqt^ 1 we have 



3(3c 3 + c 5 )Z 2 (X 2 -y 2 ) 2 
a A 



(42) 



If this value of \Q\ is 3> 1, and the Universe starts off in 
the high acceleration regime, then the whole of cosmo- 
logical history could have taken place there. In this case, 
as discussed above, some non-baryonic cosmological dark 
matter would be required. Alternatively, the decelerat- 
ing phase of cosmological expansion could push the value 
of Q into the low acceleration regime at some point, de- 
pending on parameters. One may then look to the low 
or intermediate acceleration regimes for the additional 
dust-like degrees of freedom. 

As a further possibility, we note that one could con- 
sider including extra interaction terms in the action that 
would behave like an extra dust-like contribution during 
the matter dominated era, while staying in the high ac- 
celeration regime. Indeed, such a contribution was found 
by the authors of [|[ when they considered a bimetric 
theory with the interaction term 



gig- 1 ) 



ah g a ^x. 



(43) 



One would expect that the addition of such a term in the 
present theory should have a similar effect in the | Q \ ^> 1 
regime, but it is not immediately clear what effects this 
would have on the weak-field limit. 



Low acceleration limit 



regime (with \Q\ ^> 1), it could end up in the low ac- 
celeration regime (with \Q\ <§C 1) after a suitable period 
of decelerating expansion. However, because Q depends 
on the Cj's, as well as Hubble-parameter like terms and 
I, the values of Hq and I do not necessarily tell us any- 
thing about the value of Q today. At an extreme, if the 
q's are chosen appropriately, the Universe could be in 
the low acceleration limit even before recombination, in- 
dependent of H and I. Thus, one could also envisage 
a situation where the Universe is in the low acceleration 
regime throughout its entire history. 

Now let us consider what the probes of observational 
cosmology, discussed above, can tell us about the possi- 
bility of cosmological expansion in the low acceleration 
limit. It can be seen from Table U that in this regime a 
period of cosmological expansion of the form a ~ i 1 / 2 can 
be achieved as a result of dust-like fluids being coupled 
to each of the two metrics. We therefore have the very 
unfamiliar situation of primordial nucleosynthesis being 
able to occur during a period of dust domination. How- 
ever, it can also be seen that a ~ t 1 ^ 2 can be achieved 
for a variety of other fluids coupled to the two metrics, 

. (1+4") 

including radiation, if we allow Y ~ tu+s™) . 

Now consider structure formation. It can be seen that 
a period of expansion of the form a ~ i 2 / 3 can be achieved 
if two fluids with equation of state w = —1/4 are coupled 
to each of the two metrics. Having to invent two exotic 
fluids in this way is somewhat distasteful, but again we 
have some freedom. If we are prepared to consider a sec- 
ond fluid with w = (3+16n) -1 , then we can have a ~ i 2 / 3 
being produced with a fluid of dust coupled to the first 
metric. For a theory with n = —1/8 this corresponds 
to dust and a scalar field. Various other situations can 
be read off from Equations (|3"3"]) and (|3~4"1) . One particu- 
larly interesting example is for a theory with n = —1/16. 
In this case spatial curvature dominating the first met- 
ric, and dust coupled to the second metric, also gives an 
evolution a ~ £ 2 / 3 . Here, then, the dust coupled to the 
second metric acts, in a way, as if it were cosmological 
dark matter coupled to the first, which is itself an empty, 
open universe. Clearly there is considerable scope for 
new and interesting behaviour. 



C. The intermediate regime 

One may also ask whether a dust-like contribution 
could arise from the intermediate regime, where \Q\ ~ 1. 
As only the asymptotic limits of this function are defined, 
there is considerable flexibility in terms of the function's 
transitional behaviour. For example, if we consider a 
transitional regime where M ~ \Q\ S then an analysis 
similar to that used in the low acceleration regime finds 
solutions a cx t p and X oc Y cx T q where 



As discussed above, it seems entirely plausible that 
even if the Universe starts off in the high acceleration 



2s(l - 3u>) 



(3 — 2s + 3w + 12nw — 9w — I2nw + 6sw — 9ww) 
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q = 



2s(l - 3iu) 



(3 — 2s + 3w + 12nw — 9w — 12nw + 6sw — 9ww) 



These equations reduce to (|33|) and ([34]) when s = 3/4. 
It can be seen from the above that, among very many 
other possibilities, we can arrange to have a ~ t 2 / 3 and 
Y ~ t 2 / 3 when to = w = s — 1. For a theory in which 
s = 2/3, we can therefore have the some of the effects 
of cosmological dark matter when both metrics are ef- 
fectively empty, with negative spatial curvature. This is 
intended as an example only. Without any well defined 
functional form of Ai{Q) in the intermediate regime one 
clearly has considerable freedom to achieve whatever evo- 
lution is desired. 

One final possibility arises if the entirety of cosmo- 
logical history corresponds to values of Q with a differ- 
ent sign to that which is appropriate for the weak-field 
limit. Here we have considered the functional form of 
Ai(Q) to be set by the magnitude of Q only. If one al- 
lowed its form to be different for positive or negative Q, 
then the possibility could arise that the form of A4(Q) 
that is appropriate for cosmology is not fixed at all by 
considerations of the weak-field limit. Cosmological so- 
lutions could then be considered entirely independently 
from other phenomenology. In this case one may look to 
the discussion of \Q\ ~ 1, above, for what the cosmolog- 
ical dynamics could look like for simple power-law forms 
of M{Q). 



VII. DISCUSSION 

We have studied the FRW solutions of Milgrom's class 
of bimetric theories of gravity. These theories have dif- 
ferent behaviours depending on the value of the scalar 
quantity \Q\, which is formed from the two metrics and 
various parameters of the theory. The two regimes that 
result are referred to as the 'high' and 'low acceleration' 
limits. 

We find that in the high acceleration limit the cosmo- 
logical dynamics of the two metrics essentially decouple 
from each another, and evolve in a similar fashion to the 
FRW solutions of General Relativity In this regime, the 
theory also potentially provides an explanation of the 
cosmological constant problem, as it is possible to in- 
clude a constant term in the Friedmann equations that 
is constructed from a factor of order unity divided by 
the square of the intrinsic length scale of the theory, I. 
This results in the correct order of magnitude for the cos- 
mological constant measured in our observable Universe. 
Solutions that stay in the high acceleration regime for 
their entire history, however, appear to be in conflict with 
observations, unless non-baryonic dark matter fields, or 
extra interaction terms in the action, are also included. 

In the low acceleration limit we find that simple power- 
law, matter dominated solutions also exist. There is 
considerable freedom in the form of these solutions, de- 
pending on the parameters of the theory, as well as the 



matter fields coupled to each of the two metrics. The 
solutions in this regime are, in general, quite different to 
their general relativistic counterparts, and display some 
interesting possibilities. For example, in some theories 
it is possible to have a metric that evolves like a dust 
dominated universe does in General Relativity, while it- 
self being empty and spatially open. This allows for the 
possibility of accounting for some of the effects of cosmo- 
logical dark matter via the fields coupled to the second 
metric. 
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Appendix A: Derivation of the Field Equations 

To derive the field equations, let us consider the vari- 
ations of M., a, and \f^g with respect to g ab separately. 
Firstly, varying A4 gives 



a5M = crl 2 M c 



bce f 

^L—C a bc C d ef 6g rs (44) 

+2al 2 M Q J d f 8T d ef 

where Mq = dM/dQ, and J/ f = Q ad bcef C a bc . The 
second term in the equation above can then be written 

2al 2 M Q J d ef ST d ef 
= al 2 M Q J d ef g dz (5g zf . e + Sg ez , f - 5g ef , z ) 



l z 



aM< 



(J iab) e + J (a e b )- Je a b )]j9 ab , (45) 



where = means equal up to total divergences, and where 
we have used the usual result Sg e f = —g ea gf b 8g ab . The 
quantities J ab e and J e ab , in the equation above, have had 
their indices lowered with g ab . 
We also have 



oM 5 -^ = - a -Mg ab 5g« b , 



and 



MS a 



5a 



a b 



MSg 



(ib 



(46) 



(47) 



If a = n 2n then 5a = —^ag ab 5g ab , and the last equation 
can then be written 



MSa = --ag ab MSg ab . 



(48) 



Summing (THl) . (flfj|) and (T4"8l gives the field equations 
(TT21 and shown in Section HU 
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Appendix B: The Constraint Equations 

Given the complicated form of the full second-order 
field equations, it is useful to look for the simpler, first- 
order constraint equation. To find this, consider general 
variations of the action of the form 



where 



and 



51 = / d 4 x^g~(E ab Sg ab + E ab Sg 



E a b — PG a b — S a b — T^ab — 



E a b — \ -{otGab ~ S a b ~ T^Tab) = 0. 
V 9 2 



(49) 



(50) 



If we assume that T ab and T ab contain at most first-order 
derivatives then we do not need to consider any com- 
pensatory contributions from their evolution equations 
in order to find the desired constraint equations. 

Now consider the variation of g ab and g ab due to dif- 
feomorphisms generated by the 'infinitesimal' vector field 
£ a . Firstly, let us consider the variation of Sg ab due to 
£ a , which gives 

Slug = - j ' E ab (g ac V c e+g cb V c C). (52) 

We can integrate this by parts to yield a boundary term, 
that we neglect, the Bianchi identity, V a G af> = 0, and 



(53) 



Slug = 2 / V c (£ ca )£ a . 
Now consider variations of g ab , such that 
SIuS = J E ab {^V c g ab - g ac V c ( b - g cb V c n (54) 
= / (UV c g ab )Eab + HaVc(E b a g bc + E\g cb )). 



Collecting terms then gives 

5 h = h,Sg + J Ug ( 55 ) 
= / U(V a g cb )E cb + W c (2E ca + E b a g bc + E\g cb )). 

I 



Finally, let us define the tensor J ca = 2E ca + E b a g bc + 
E a b g cb . If J 0a had second-order derivatives, then V c J ca 
would contain third-order derivatives. However, we know 
from the vanishing of the variation of the action under 
diffeomorphisms, 5I ( = 0, that V c J ca = -(V a g cb )E cb . 
Now, due to the structure of the theory we can see 
that E cb contains at most second-order derivatives, so 
the right-hand side of this equation can contain at most 
second-order derivatives. We can therefore conclude that 
J 0a contains up to first-order derivatives only. 

For a given foliation of hypersurfaces with normal n a , 
and metric h ab = g ab + n a n b , we then find that the con- 
straint equations are given by 



(51) and 



J ab n a h b r = 0. 



J ab n a n b = 



(56) 



(57) 



The first of these is trivially satisfied by FRW geometry, 
where the hypersurfaces are taken to be surfaces of con- 
stant t, while the second gives the constraint equations 
displayed in Appendix C, and Section HVl 



Appendix C: The Field Equations with FRW 
Geometry 

Here we will write explicit expressions for some of the 
quantities that appear in the field equations (1121) and 
(1T31) . for the FRW geometries CE?]) and (H]). Firstly, we 
can immediately write 



\/l — kr 2 



VI - kr 2 XY 3 



(58) 



We can also write a relatively simple expression for Q. 
To do this we first define five new quantities via 



Q = ctQW + c 2 QW + c 3 Q^> + c 4 QW + ^QW (59) 

such that Q« = l 2 5 f J d b g ce C a cd C b ef , with the other 
four defined mutatis mutandis. We can then write 



(1) 


I 2 

a 2 


4- 

a z 


aX 
f \l 


X 2 
X 2 


Y 2 


(2) 


I 2 
a 2 


' a 2 


_ 2 aX 
aX 


aY 


X 2 
X 2 



(k - kfr 2 



(1 - fcr 2 )(l - kr 2 ) 2 
XY (fc-fc) 2 r 2 (3 



2((1 -kr 2 ) + 2(1 - kr 2 )) YY ( Y a 

'X 2 \ Y~a 



(1 - kr 2 ) 



2kr 2 



XY 



{1 - kr 2 ){\ - kr 2 f 
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({1 - kr 2 ) + 2(1 - kr 2 )) YY X Y 



(1 - kr 2 ) 



X 2 \ X 



+ 3— -4- 



Y 



)(3) 



,-,2 



2 \2 



a 2 aX X 2 (k-k) 2 r 2 (Z-2kr 2 ) 
4 4 1_ 

^ 2 al X 2 



a' 



(l-fcr 2 )(l-/cr 2 ) 2 
((1 - /cr 2 ) + 2(1 - fcr 2 )) rr 



(1 - kr 2 ) 



X 2 



a n X ((1 - kr 2 ) + 2{l - kr 2 )) YY 

4 v ~- X 2 



Q (4) 



16- 



a X (1-fcr 2 ) 
(k - kfr 2 



aX aY X 2 XY Y 2 
5-— + 24-— - — - 6 — - - 9— + - 

aX aY X 2 XY Y 2 (I - kr 2 ){l - kr 2 ) 



aX 



aY X< 



Y 2 



lO^r + 2-— + 12-- - —7 - 6— r 



(/c-fc) 2 r 2 (l + 2(l-/cr 2 ) 2 ) 



a y x 2 y 2 (i - fcr 2 )(i - ifer 2 ) 2 

((1 - /cr 2 ) + 2(1 - kr 2 )) YYa ((1 - kr 2 ) 2 + 2(1 - kr 2 ) 2 ) Y 2 Y 2 



(1 - fcr 2 ) 



X 2 a 



(1 - fcr 2 ) 2 



X 2 



When k — k (motivated in Section |V| , the constraint equation can then be written for arbitrary c\ , c 2 , C3 , c 4 and C5 
in relatively simple form: 



2crA^ c 



r 3 /y 2 



x 2 



a 2 x \ y 2 +fc y 2 



a 2 7Wa 87rG/9a 2 87rGpXy 3 



3/ 2 



3a 2 



2(ci - 2c 2 - 2c 3 - 8c 4 - 5c 5 )% + (2ci + 5c 2 - 4c 3 + 8c 4 + 2c 5 )-^ 



X 2 3((c 2 + 8c 4 + 4c5)A- 2 + 2(-ci + 2c 3 + c 5 )y 2 ) a y 

x 2 + x 2 a y 

3((c2 + 2c 4 )A" 2 - 2c 3 y 2 ) X y 3((c x + 3c 4 + 2c 5 )A" 4 - 2 Cl X 2 y 2 + (3c 3 + c 5 )y 4 ) y 2 



-(ci + c 2 + c 3 + c 4 + c 5 )-^ 



A" 2 
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A" 4 



y 2 



where cr and Q are given by (|58|) and J59J), above. The special case of the 'concrete simple theory' of Milgrom is 
given in Section IIV1 Due to their length, we choose not to display the evolution equations here, which make up the 
remaining part of the field equations. 
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